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Using a multistep renormalization group method, we study the low-temperature phases of the 
interacting one- dimensional (ID) electron gas coupled to phonons. We obtain analytic expressions 
for the weak-coupling quantum phase boundaries of the ID extended Holstein-Hubbard model and 
the ID extended Peierls-Hubbard model for general band-filling and phonon frequency. Away from 
half-filling, the phase diagrams are characterized by a delicate competition between spin density 
wave, charge density wave, and superconducting orders. We study the dependence of the ground 
state on the electron-phonon (el-ph) and electron-electron (el-el) coupling strengths, the screening 
length, electron bandwidth, phonon frequency, doping, and type of phonon. Unlike the case in Fermi 
liquids, in ID the el-ph coupling is strongly renormalized, often to stronger values. Even when the 
bare phonon-induced attraction is weak compared to the bare el-el repulsion, a small amount of 
retardation can cause the renormalized el-ph interaction to dominate the problem. We find cases 
in which a repulsive el-el interaction enhances the superconducting susceptibility in the presence 
of a retarded el-ph interaction. The spin gap and superconducting susceptibility are found to be 
strongly dependent on the deviation from half-filling (doping). In some cases, the superconducting 
susceptibility varies nonmonotonically with doping and exhibits a maximum at a particular doping. 
For a quasi-lD array of weakly coupled, fluctuating ID chains, the superconducting transition 
temperature T c also exhibits a maximum as a function of doping. The effect of changing the ion 
mass (isotope effect) on T c is found to be largest near half-filling and to decrease rapidly with doping. 
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I. INTRODUCTION 

Recent experiments 1,2 in the high-temperature super- 
conductors that suggest a strong and ubiquitous electron- 
phonon (el-ph) interaction have lead to an increased in- 
terest in strongly correlated el-ph systems. Of particular 
interest is the interplay between the repulsive, instanta- 
neous Coulomb interaction and the attractive, retarded 
interaction mediated by phonons, in the context of var- 
ious phase transitions (instabilities) of the system. This 
interplay is quite simple and well understood in a two- 
or three-dimensional Fermi liquid, and serves as a foun- 
dation for the BCS theory of superconductivity. How- 
ever, the high-tcmperature superconductors, as well as 
some other materials of current interest, exhibit mani- 
festly non-Fermi liquid behavior 

Compared to conventional metals, much less is known 
about the influence of an el-ph interaction in non-Fermi 
liquids. One such system, where much analytic progress 
can be madej&^iiMiiiSi is the interacting one-dimensional 
electron gas (1DEG), where the system is strongly corre- 
lated even for weak interactions. In contrast to higher 
dimensions, in one-dimension the el-ph interaction is 
strongly renormalized. In other words, the effective el- 
ph interaction at low energies is not simply given by the 
bare, microscopic coupling constants, which makes the 
physics both very different and very rich. Furthermore, 
the renormalizations of the el-ph interactions are affected 
by direct electron-electron (el-el) interactions. It is un- 
known how much of the physics of the 1DEG is generic to 
other non-Fermi liquids; but one may hope that certain 
features are generic, in which case the 1DEG may serve 
as a paradigmatic model for a broader class of systems. 



In addition, there are many reasons to study the 1DEG 
coupled to phonons for its own sake, including the im- 
portance of el-ph interactions in the quasi-lD conducting 
polymers^ii4 and the possibility that they play an im- 
portant role in the quasi-lD organic conductors^ It may 
also be the case that holes in the high-temperature super- 
conductors live in quasi-lD due to stripe correlations. 16 
In any case, since in ID a small bare el-ph coupling can 
be renormalized to substantially larger values, the el-ph 
interaction is important to consider. 

In the present paper, we use a multistep renormal- 
ization group (RG) procedure to comprehensively study 
the zero-temperature phase diagram of the spinful 1DEG 
coupled to phonons, treating the el-el and el-ph interac- 
tions on equal footing. The same technique is employed 
to compute the doping dependent superconducting sus- 
ceptibilities, charge density wave susceptibility, and iso- 
tope effects. In a separate paper ^ we have studied the in- 
fluence of the el-ph interaction on the electron dynamics 
of an interacting 1DEG, expressed via the single particle 
spectral function. The strategy in the present paper is 
to start with a microscopic electron-phonon model with 
many parameters (el-el interactions, electron bandwidth, 
el-ph interactions, and phonon frequency), and then to 
integrate out high-energy degrees of freedom to produce 
a low-energy effective field theory with a known phase 
diagram-the continuum IDEG-whose only parameters 
are the renormalized el-el interactions and bandwidth. 

If the system is far from commensurate filling, we em- 
ploy the "two-step RG" procedureii2ii£ii£ln this method, 
one set of RG equations governs the flow of the coupling- 
constants for energies u> greater than the phonon fre- 
quency wq, while a second set governs the flow for loq > u>. 
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If, as usual, the Fermi energy Ep > ui Q , the first step is 
to integrate out degrees of freedom from Ep to u>q using 
the microscopic coupling constants as initial values. The 
resulting renormalized couplings are then used as initial 
values in the second stage of RG flows, to integrate out 
degrees of freedom from luq to some low-energy scale. 
We also use the two-step RG technique to study systems 
at half-filling (the only difference being the inclusion of 
Umklapp scattering). Near, but not equal to half- filling, 
there are three steps to the RG transformation. Defin- 
ing /i as the chemical potential relative to its value at 
half-filling, for \i > uj n these steps are Ep > uj > fj,, then 
/j > uj > ujq, and finally ojq > ui. This "three-step RG" 
technique allows us to study the continuous evolution of 
the phase diagram with doping. 

We study the continuum limit of two microscopic 
models of interacting, spinful ID electrons coupled to 
phonons: the extended Holstein-Hubbard modeli^ and 
extended Peierls-Hubbard model. While others have em- 
ployed RG techniques to study the 1DEG coupled to 
phonons ji2*ii*i2iiLi& and indeed some of the qualitative 
results of the present paper have been known for some 
time, the aforementioned models have not been explored 
in detail with multistep RG. Numerical calculations of 
their phase diagrams have been mostly limited to sim- 
plified models that contain spinless electrons, infinite ion 
mass, or zero el-el interactions ,SL2i These studies have 
also been mostly limited to specific values of the band 
filling, especially half- filling. 22 

We obtain analytic expressions for the low- 
temperature phase boundaries of the aforementioned 
models. Since the RG procedure is perturbative (one- 
loop), our results are only accurate for interactions that 
are small compared to the bandwidth, but are valid for 
any relative strength of the el-ph and el-el interactions. 
Corrections at stronger couplings are expected to be 
smooth and should not make large qualitative changes 
to the results, as long as the interactions are not too 
large. The method properly takes into account the 
quantum phonon dynamics, and is therefore used to 
study phonons of nonzero frequency. 

One question we address is whether superconductiv- 
ity can exist in realistic quasi- ID systems in which the 
bare el-el repulsion is stronger than the bare attractive 
interaction mediated by the el-ph coupling. It is known 
that in the absence of el-ph interactions, a repulsive el- 
el interaction in a single-chain 1DEG is always harmful 
to superconductivity, as is the case in higher dimensions. 
However, we find that for the 1DEG coupled to phonons, 
increasing the el-el repulsion can in some specific cases 
enhance superconductivity. Moreover, with even a small 
amount of retardation present, it is possible for a 1DEG 
to have a divergent superconducting susceptibility even 
when the bare el-el repulsive is much stronger than the 
bare el-ph attraction. (In three dimensions, this is only 
possible with substantial retardation.) 

In ordinary metals, the observation of an isotope ef- 
fect on T c was crucial in the development of the BCS 



theory that describes the superconductivity of Fermi liq- 
uids. The isotope effect exponents ap c = —d\nT c /dlnM 
and a a = — <i In A/dln M have the universal value of 1 /2 , 
where M is the ion mass and A is the superconducting 
gap. 

In contrast, in the cuprate high-temperature supercon- 
ductors, both T c and ap c are strongly doping-dependent. 
Despite the fact that the superconducting gap is a mono- 
tonically decreasing function of increasing doping, T c 
varies nonmonotonically with doping, exhibiting a maxi- 
mum at "optimal doping." For dopings well below opti- 
mal, the isotope effect on T c is quite large: ap a ~ 1. As 
the doping increases, ap c decreases, usually dropping be- 
low 0.1 near optimal doping^ The isotope effect on the 
so-called pseudogap has the opposite sign as the isotope 
effect on T C H The origin of these highly unconventional 
isotope effects remains one of the many unsolved myster- 
ies of high-temperature superconductivity. 

Since the 1DEG is perhaps the only presently solvable 
non-Fermi liquid, it is worth computing isotope effects in 
this system to try to gain insights on isotope effects in 
unconventional superconductors. We compute ap c for a 
quasi-lDEG coupled to phonons, under the assumption 
that charge density wave order is dephased by spatial 
or dynamic fluctuations of the ID chains. 25-26 For most 
choices of the parameters, ap c is larger than the BCS 
value at small dopings, then drops below 1/2 as the dop- 
ing is increased. We show that the quasi-lDEG coupled 
to phonons displays a strongly doping-dependent T c that 
can exhibit a maximum as a function of doping. This 
behavior occurs despite the fact that the pairing energy, 
determined by the spin gap A s , is a monotonically de- 
creasing function of increasing doping. We also compute 
the isotope exponent «a 3 = — rflnA s /c?lnM and find 
cka s < 0, which in most cases is the opposite sign as q^. 

The rest of this paper is organized as follows. SectionlTTI 
defines the microscopic models. Section ITTT1 presents our 
results for the phase diagrams, without derivation. In 
Section lTvl we give our results for the doping dependence 
of the superconducting susceptibility and isotope effects, 
again without derivation. In Section we compare our 
analytical results for the phase diagrams to numerical 
work of other authors. Scction lVll discusses the RG flows 
of the coupling constants and contains a mathematical 
derivation of all the results in the previous sections. In 
Section IVIII we summarize the results qualitatively and 
make some concluding remarks. 



II. MODELS OF ID ELECTRON-PHONON 
SYSTEMS 

The ID extended Peierls-Hubbard (Pei-Hub) model is 
defined by the Hamiltonian 

Wpei-Hub = ~<y^[l - + Q»)](4.a C »+l,tr + H.C.) 

Z,<T 



where the el-el interaction part is that of the extended 
Hubbard model: 
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Here, c| CT creates an electron of spin er on site i, a- cre- 
ates a phonon of frequency lvq between sites i and i + 1, 

rij = X)(T n i,o" = EtT c l cr c i,cr, and 7 is the dimensionless 
el-ph coupling constant. This model, in the absence of 
Tiuv, is an approximation to the model of Su, Schrieffer, 
and Heeger— (SSH). Including extended Hubbard inter- 
actions, the SSH model is 

^SSH-Hub — 

~ E ^ ~ 7^+1 ~ ( C lcr C i+h<? + H - C -) 



where vp = 2tsin(kp), V = — V cos(2kp), Ep > iv is 
a high-energy cutoff for the RG theory on the order of 
the Fermi energy, and 2kp/ir is the average number of 
fermions per site. The el-ph coupling parameters Ap C i 
and AhoI are defined such that, in the absence of el-el 
interactions, the spin gap is given by A s oc exp(— a/A), 
where a = 1/2 for half-filling, a = 1 for incommensurate 
fillings, and A stands for Ap c i or AhoI, depending on the 
model. (In Section IVII we give the result for A s is the 
presence of el-el interactions.) The method we employ 
yields phase diagrams that are accurate for A, U, V <C 1. 
We expect that the technique is qualitatively accurate 
when these couplings are of order 1. It should be clear 
that, in the present paper, we have set the lattice param- 
eter and Planck's constant equal to 1. 
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Here, acoustic phonons with spring constant k couple to 
electrons by modifying the bare hopping matrix element 
t by the el-ph coupling strength 7 times the relative dis- 
placements Ui+i — Ui of two neighboring ions of mass M. 
If we approximate the acoustic phonon as an Einstein 
phonon of frequency 2^/k/M = ujo, the SSH-Hub model 
reduces to the Pei-Hub model. This is a good approx- 
imation since, in the SSH model, the el-ph interaction 
vanishes at zero momentum transfer. The el-ph coupling 
constants of the models are related via 7 = 7 t\/2MuiQ. 

The ID extended Holstein-Hubbard (Hoi-Hub) model 
is defined by the Hamiltonian 



T^Hol-Hub = 

i,a i 

+ g\/2Moj 2J Qi^i + Hi 
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In this model, a dispersionless optical phonon mode with 
vibrational coordinate qi and frequency ojq couples to the 
local electron density with el-ph coupling strength g. Un- 
like the Pei-Hub model, this model contains equal parts 
backward scattering (momentum transfer near 2kp) and 
forward scattering (momentum transfer near 0) el-ph in- 
teractions. Another difference is that the el-ph interac- 
tion is site centered (diagonal) in the Hoi-Hub model ver- 
sus bond centered (off-diagonal) in the Pei-Hub model. 
In the present paper, we only explicitly discuss the case 
of repulsive el-el interactions (U, V > 0); however, the 
mathematics remains valid for attractive ones. 

For convenience we define the following dimensionless 
quantities: 
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III. RESULTS FOR THE PHASE DIAGRAMS 

In this section, we present our main results for the 
phase diagrams, without derivation. More discussion of 
the method and a detailed derivation are given in Section 
IVII where we obtain explicit expressions for the phase 
boundaries. 



A. Incommensurate filling 

Below we present zero-temperature phase diagrams in 
the incommensurate limit, which corresponds to fi ~ Ep. 



1. Transition to the spin-gapped phase 

For the ID extended Hubbard model without el-ph 
interactions, the low-energy properties are described by 
the spin-charge separated Luttinger liquid (LL) as long 
as V < U/2. The low energy properties of this gap- 
less, quantum-critical state of matter can be described 
by a bosonic free field theoryi^Si Since the quasiparti- 
cle residue vanishes, the LL is by definition a non-Fermi 
liquid; there are no elementary excitations with the quan- 
tum number of an electron (or a hole). 

In the absence of el-ph interactions, a spectral gap de- 
velops in the spin sector if V > U/2, which leads to quite 
different physical properties than the LL. This non-Fermi 
liquid phase is termed a Luther-Emery liquid^ (LEL). 
However, in nature one typically expects V < U/2, so 
that some additional physics is needed to created a spin 
gap. For incommensurate fillings, the charge sector is 
gapless. 

We now study the effects of the el-ph interaction on the 
phase boundary between the LL and spin-gapped LEL. 
In Fig. ^ we show a phase diagram in the V-U plane, for 
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FIG. 1: Phase diagram in the V-U plane showing the 
A dependence of the phase boundary between the gapless 
Luttinger liquid (LL) phase and the spin-gapped Luther- 
Emery liquid (LEL) phase, for an incommensurate 1DEG 
with Ef/(±>o = 5 (dashed lines) and Ef/ujq = 10 (thin solid 
lines). The thick line shows the phase boundary at A = for 
any Ef/loo- A stands for either A Ho i or Ap c i, depending on 
the model. 
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FIG. 2: Phase diagram in the XIq-UIq plane for an incommen- 
surate 1DEG, showing the dependence on V/U of the phase 
boundary separating the LL and LEL phases. A stands for 
either Ah i or Ap c i, depending on the model. 
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FIG. 3: Phase diagrams of the incommensurate extended 
Hoi-Hub model with (a) V = and (b) V/U = 0.2, and of 
the incommensurate extended Pei-Hub model for (c) V = 
and (d) V/U = 0.2. For all diagrams, E F /u = 10. To 
the right of the thick line, the system is spin-gapped. The 
most divergent susceptibility is shown without parenthesis, 
while parenthesis indicate a susceptibility that diverges less 
strongly. SDW stands for 2k F spin density wave, CDW stands 
for 2k,F charge density wave, SS stands for singlet supercon- 
ductivity, and 4kF stands for 4&f charge density wave. 



2. Competition between spin density wave, charge density 
wave, and superconducting instabilities 



various fixed values of A and Ef/loq (this phase bound- 
ary is the same for the Pei-Hub and Hoi-Hub models). 
We see that a retarded el-ph interaction dramatically in- 
creases the stability of the LEL phase relative to the LL 
phase. The phase boundary is very sensitive to the re- 
tardation parameter Ef/loq, with higher values favoring 
the LEL phase. For the case of an unretarded el-ph inter- 
action (luq > Ep), which is not typical in real materials, 
a spin gap can only occur for A > U — 2V. However, 
just a modest amount of retardation makes a spin gap 
possible, even, in some cases, for A <C U — 2V. As we 
show in Section IVII this is due to the renormalization of 
the backscattering el-ph interaction toward stronger val- 
ues. The figure shows that poorly screened interactions 
(high V/U) favor the LEL. The dependence of the phase 
boundary on V/U is shown explicitly in Fig. [21 which 
presents a phase diagram in the Xlg-Ulo plane. This di- 
agram shows that scaling eventually carries one to the 
case in which an infinitesimal A causes a spin gap (in 
other words, for Ep 3> too, the system is spin-gapped for 
infinitesimal A). 



Below, we explore the many ordering instabilities 
present in the system. 2kp spin density wave order 
(SDW), 2kp charge density wave order (CDW), and sin- 
glet superconductivity (SS) can all compete at zero tem- 
perature. A divergent charge density wave susceptibil- 
ity with Akp periodicity (labeled in phase diagrams as 
tl Akp") is also possible. It is important to note that 
since long range order is forbidden in an incommensurate 
ID system, the phase diagrams below actually consist 
of identifying instabilities with divergent response func- 
tions. However, for a quasi-lD array of weakly coupled 
chains, interchain coupling allows for true broken sym- 
metry order at low temperature. 

In Fig. |3 we present phase diagrams in the £/-Aho1 
and U-Xp Q i planes, for Ep/luq — 10. In these diagrams, 
we show phase boundaries between regions where vari- 
ous ordering fluctuations have divergent susceptibilities 
in the low-temperature limit. The susceptibility that di- 
verges most strongly i.e., dominates, is shown without 
parenthesis. If a second susceptibility diverges, but less 
strongly, it is termed "subdominant," and is shown in 
parenthesis. The thick solid line is the LL-LEL transi- 
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FIG. 4: Phase diagrams of the incommensurate extended 
Hoi-Hub model, (a) and (b) Ah i = 0.2; (c) and (d) AhoI = 
0.3. (a) and (c) V/U = 0; (b) and (d) V/U = 0.1. Parenthesis 
indicate a subdominant susceptibility. The region with dom- 
inant SDW is the LL phase, while the rest of the parameter 
space is a LEL. 
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tion line; the LL phase is present to left of this line and 
the LEL phase to the right. 

For repulsive el-el interactions, the entire LL phase, for 
either model, is dominated by SDW fluctuations, with a 
slightly weaker CDW susceptibility. The LEL phase is 
more complex. For the extended Hoi-Hub model, dom- 
inant SS order is possible provided the el-el repulsion is 
weak enough and AhoI is neither too weak nor too strong. 
For the Pei-Hub model with repulsive el-el interactions, 
a phase with dominant SS is impossible due to the ab- 
sence of el-ph forward scattering. Therefore, generally 
speaking, an optical phonon is more favorable to super- 
conductivity than an acoustic one. In both models, there 
is a large region, for intermediate values of A, with domi- 
nant CDW and subdominant SS. At high values of A, SS 
is no longer divergent; in this region a 2k f CDW domi- 
nates and a Akp CDW is subdominant. We must point 
out that the dashed line is not expected to be quanti- 
tatively accurate, since the method is a weak-coupling 
one. Note that the phase with dominant SS is strongly 
suppressed by poor screening (large V/U). 

In Fig. 21 we study the dependence of the ground state 
on Ef/wo, for the Hoi-Hub model, by showing phase dia- 
grams in the Ef/ojq-U plane. We see that high values of 
Ep/u create a spin-gapped phase with dominant CDW. 
Low values of Ep/too create a SS-dominated LEL for low 
U, and a LL for high U. In other words, an electron 
bandwidth that is too big is harmful to superconductiv- 
ity! Note that for moderate Ep/uiQ, the system lies in the 



region with dominant CDW and subdominant SS, which 
extends from U — to quite large values of U. Therefore, 
even when the bare interactions are predominantly repul- 
sive (U ^> Ahoi), it is still possible for the system to have 
a divergent superconducting correlation. It is possible 
for this to occur even with a small degree of retardation 
like Ep/uJo ~ 5 [see Fig. Eld)]. The phase diagram of 
the Pei-Hub model is similar to Fig. 0] except that the 
phase with dominant SS is removed, and the dashed line 
is shifted to slightly lower Ep/ujQ. 

It is worth pointing out the intriguing possibility that, 
for a quasi-lD system with dynamically fluctuating ID 
chains, or even for chains that exhibit transverse spatial 
fluctuations, CDW order is easily dephased, while the su- 
perconducting instability is noti2& If this is the case, then 
it is possible for the system to support superconductiv- 
ity even for the physically realistic case of U ^> AhoI, and 
without the large amount of retardation that is required 
in 3D. 

A number of authors have computed the phase dia- 
gram of the spinless Holstein model at half-filling, in 
the absence of el-el interactions (see Fig. To fa- 

cilitate comparison between this model and the spinful 
incommensurate extended Holstcin-Hubbard model, Fig. 
shows our phase diagram for the latter model, in units 
similar to Fig. El 
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3. Enhancement of superconductivity by repulsive 
interactions 

The intuitive notion that repulsive interactions sup- 
press superconductivity at weak-coupling, while always 
true in a Fermi liquid, does not always hold for the 
1DEG coupled to phonons. In the 1DEG, the poten- 
tially strongly divergent part of the singlet supercon- 
ducting susceptibility at temperature T <C u>o is Xss oc 
AsT 1 ^ ~ 2 . Interactions in the spin channel determine 
A s , while charge interactions renormalize the effective 
Luttinger parameter K° s away from its noninteracting 
value of 1 (see Section IVl) . Clearly, SS is enhanced by 
an increase in A s or an increase in K c . In many cases, 
increasing the el-el repulsion causes both A s and K c to 
decrease. However, below we discuss cases in which one 
of the two parameters is increased by el-el repulsion; de- 
pending on how much the second parameter is reduced, 
Xss may be enhanced. 

In the absence of el-ph interactions, an increase in U 
always decreases K c and A s , and therefore suppresses 
superconductivity. In a ID el-ph system, this is often 
the case as well; however, there are also cases in which 
increasing U can cause K° s to be increased. (The tech- 
nical reason is a decrease in the effective el-ph backscat- 
tcring interaction-see Section WH ) Then, if A s is not 
reduced too much by the increase in U, it is possible for 
SS to be enhanced. An example of this can be seen in 
Fig. EJa) orBJc). If we start at A = 0.35 and [7 = 0, 
and increase U while holding A fixed, we cross from a 
region without divergent SS to a region with divergent 
SS. This phenomenon can also be seen in Figs. EI a ) and 
0Jc) if the right value of Ef/ojq is chosen [such as, for 
example, Ef/ujq = 200 in Fig. Ufa)]. However, if we 
set V/U > 1/6 and hold this ratio fixed while increasing 
U, superconductivity is never enhanced at weak coupling 
[see Figs. Etb),EId), and Section IVTj. 

In Fig. H3 we further investigate the enhancement of 
SS by repulsive interactions by plotting the dimensionless 
superconducting susceptibility 
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XSS = TTfFXSS 



(A S /E F )(T/E F )V K ?- 2 (6) 



versus U, at fixed A = 0.3, V = 0, and T/uj = 0.01. For 
small U, xss increases with increasing U. However, as U 
is increased further, K^ s stops increasing as rapidly, and 
the decreasing A s causes Xss to drop back down. 

It is also possible to enhance superconductivity, in 
some cases, by increasing the nearest-neighbor Hubbard 
repulsion V, while holding U fixed. This causes a renor- 
malization of the el-ph backscattering toward stronger 
coupling, which results in an increase of A s and a de- 
crease of K% . Depending on the competition between 
these two effects, xss may (or may not) be enhanced. An 
example of a case in which xss is enhanced can be seen in 
Figs. OJa) and^Jc). There, if one begins at A = 0.1 and 
U = 0.2, then increases V/U from to 0.2 while hold- 
ing A and U fixed, the system moves from a LL phase 
without divergent SS, to a LEL phase with divergent SS. 
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FIG. 6: Dependence of the superconducting susceptibility 
Xss (thick line) on the Hubbard repulsion U, for an incom- 
mensurate system with V = 0, fo/uo = 10, and T/ujo = 0.01. 
(a) is for the Hoi-Hub model with AhoI = 0.3, and (b) is for 
the Pei-Hub model with Ap c i = 0.3. The thin line is A s /ujo 
and the dashed line is the effective Luttinger exponent K^ S ■ 
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FIG. 7: Dependence of xss (thick line) on the nearest- 
neighbor repulsion V, for an incommensurate system with 
U = 0.3, E f /loo = 10, and T/wo = 0.01. (a) is for the Hol- 
Hub model with AhoI = 0.2, and (b) is for the Pei-Hub model 
with Ap c i = 0.2. The thin line shows A s /cjo and the dashed 
line is the effective Luttinger exponent K° s . 



In Fig. we explore this phenomenon further by 
plotting xss versus V at fixed U — 0.3, A = 0.2, and 
T/uio = 0.01. For small V, xss is enhanced by increas- 
ing V, due to the increase in A s . However, as V is in- 
creased further, eventually the rapidly decreasing K° 
begins to overwhelm the effect of the increasing A s , and 
Xss decreases. In this case, optimizing superconductivity 
therefore requires a fine tuning of V/U. 

We have pointed out exceptions to "rule" that el-el re- 
pulsion suppresses superconductivity at weak coupling, 
in order to illustrate, as a point of principle, the dramati- 
cally different physics that governs ID el-ph systems com- 
pared with a Fermi liquid coupled to phonons. It is worth 
briefly discussing some prior works on this topic. Al- 
though the RG flow equations in Ref. [I(J are correct, the 
authors implied that an enhancement of superconductiv- 
ity by repulsive interactions is a generic feature of the 
1DEG coupled to phonons, while Ref. ^3 concluded that 
repulsion always suppresses superconductivity. Our re- 
sults indicate that both works overstated things; indeed, 
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FIG. 8: Dependence of the transition line between spin- 
gapped and non-spin-gapped phases on the doping parameter 
8 and on the electron-phonon model, for V = 0. For n < loq 
(which includes the half-filled case n — 0), the transition line 
is denoted by 8 — 0. The incommensurate limit (fi ~ Ef) is 
labeled by 5 = 1. 8 = 1/3 and 2/3 are intermediate dopings, 
with (i/ujo — (Ef/ujo) 5 - Here, A stands for Ap c i in the Pei- 
Hub model (dashed lines) and for AhoI in the Hoi-Hub model 
(dashed-dotted lines). For 8=1 (solid line), the transition 
line is the same for either model. 
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FIG. 9: Doping dependence of the Hoi-Hub phase diagram 
for V = and Ef /loo = 5. The value of 8 is given above each 
plot. 



both situations are possible, depending on the choice of 
parameters. The disagreement between Ref. ^3 an d Ref. 
Il^ was caused, in large part, by the fact that Ref. fo- 
cused purely on the effect of the el-ph interaction on A s , 
while Ref. [l3 focused purely on the effect of the el-ph 
interaction on K° s . Above, we have correctly taken into 
account that the el-ph interaction affects both A s and 
K° s , both of which in turn affect xss- 



B. Near half-filling 

If we reduce the doping level of an incommensurate sys- 
tem (i.e., move closer to half- filling), the LL-LEL phase 
boundary is influenced in opposite ways for the Hoi- Hub 
compared to the Pei-Hub model. This is because the on- 
site Hubbard repulsion is in direct competition with the 
attractive on-site el-ph interaction of the Holstcin model, 
while it cooperates with the attractive bond centered in- 
teraction of the Pei model. Therefore, the spin gap is en- 
hanced by proximity to half-filling for the Pei-Hub model, 
and reduced for the Hoi-Hub model. 

We illustrate this in Fig. |Hl which presents a phase 
diagram in the UIq-^Iq plane, for V — 0. This diagram 
shows the dependence of the LL-LEL phase boundary on 
the doping parameter 



ka.(E F /wo)' 



(7) 



Assuming the charge gap A c < fi, the actual doping con- 
centration x relative to half-filling is related to S by 
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FIG. 10: Doping dependence of the Pei-Hub phase diagram 
for V = and Ef/ujq = 5. The value of 8 is given above each 
plot. 



where v c w vf is the charge velocity [see Eq. i[41|l]. The 
incommensurate limit 8 = 1 is shown previously in Fig. 

and Fig. Ofb) as a thick solid line. As we move 
closer to half-filling by lowering 6, the LL-LEL transition 
line for the Pei-Hub model (dashed line) moves toward 
lower values of A, while the transition line for the Hoi- Hub 
model (dash-dotted line) moves toward higher values. In 
the weak-coupling limit assumed here, for fi < lu , the 
LL-LEL transition line is independent of [i. Therefore, 
this phase boundary is the same for half-filling (/i = x = 
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FIG. 11: Phase diagrams for Ef/u)q — 5 in the X-8 plane. 
Plots (a) and (b) are for the Hoi-Hub model; (c) and (d) are 
for the Pei-Hub model. Plots (a) and (c) are for U = 0.1; (b) 
and (d) are for U = 0.4. For all plots, V = 0. 



0) as for (5 = 0. 

The doping dependencies of the other phase bound- 
aries are shown in Figs. 151 and ITUl for the range u>o < 
fi < Ep, Figure[3]illustrates that for the Hoi- Hub model, 
proximity to half-filling strongly suppresses the phase 
with dominant SS. For both models, moving toward half- 
filling increases the stability of the phase with subdom- 
inant AUf charge density wave, at the expense of the 
phase with subdominant SS, especially for the Pei-Hub 
model (Fig. [TDJ). Note that, especially in Fig. EH SS is 
divergent for a range of parameters such that U 3> Ap i, 
despite the low value of Ef/ojq = 5. Figure ITTI shows a 
different slice of the phase diagram by showing plots in 
the A-c5 plane. 



C. Half-filling 

At half-filling, for repulsive el-el interactions, a charge 
gap is present for the entire phase diagram. Divergent 
SS is eliminated at half-filling. The region without a 
spin gap is a SDW, while the spin-gapped region is an 
ordered CDW. We show the half-filled phase diagram, 
for both models, in Fig. ^1 for several values of Ef/ojq 
and V = 0. The difference between the two models is 
substantial, with the Pei-Hub model favoring CDW more 
than the Hoi-Hub model. In Fig. I12f b). we also draw a 
dashed line defined by Ap c ; = U, which is the transition 
line predicted by mean-field theory. Such a treatment is 
known^S to be quite inaccurate for the Pei-Hub model, 
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FIG. 12: Phase diagram of the (a) half-filled Hoi-Hub model 
and (b) half-filled Pei-Hub model, both with V = 0. The 
dashed lines are for Ef/ujq = 5 and the solid lines for 
Ef /u)o = 10. For both diagrams, a spin gap is present for the 
CDW phase, while a charge gap is present everywhere in the 
phase diagram. In (b) the dotted line, defined by Ap e i = U, is 
the inaccurate result for the phase boundary from mean-field 
theory. 



as demonstrated in the figure, since it does not take into 
account the dramatic renormalization of the backscatter- 
ing el-ph interaction to stronger couplings. Note that the 
Pei-Hub phase diagram is very sensitive to Ep/uio, with 
high values favoring a spin-gapped CDW, while the Hol- 
Hub phase diagram is only weakly dependent on Ef/ujq. 
It is interesting that in the Pei-Hub model [Fig. I12f b)]. 
there is a maximum value of the critical el-ph coupling of 
about 0.149/^o, which occurs at U ~ 0.411/Zo- (In other 
words, for Ap c i > 0.149/Zn, the system is an ordered CDW 
for any U.) 



IV. DOPING DEPENDENCE OF THE 
SUPERCONDUCTING SUSCEPTIBILITY AND 
ISOTOPE EFFECTS 

In this section, we study the strong doping dependen- 
cies of the spin gap, superconducting susceptibility, CDW 
susceptibility, and isotope effects. 

Examining the phase diagram in Fig. IllT d). we can de- 
duce an interesting nonmonotonic dependence of the SS 
susceptibility on 5. For moderate values of Ap e i, for ex- 
ample, near Ap e i « 0.2, \SS is not divergent near 8 = 0, 
where only the 2Uf and Akp CDW susceptibilities di- 
verge, nor is it divergent near 6 = 1, where the system 
is in the gapless LL phase. However, \SS is divergent 
for a certain range of moderate 5. Therefore, in these 
cases, at fixed T <C A s , xss must exhibit a maximum as 
a function of d at some intermediate value of 5. 

This maximum in xsSj which occurs in both models, 
is shown explicitly in Figs. H3l and lT4l where we plot xss 
(thick solid line) versus 6 at T/u>o = 0.1, for representa- 
tive parameters. The cause of the nonmonotonic doping 
dependence is the different doping dependencies of A s 
and K^ B . A s decreases with increasing doping, which 
acts to reduce xsSj while K^ s increases with increasing 
doping, which acts to increase xss- These two effects 



9 



"compete" with each other and can cause a maximum 
at some "optimal" value of the doping that depends on 
the interaction strengths. The dimensionless 2kp CDW 
susceptibility 

Xcdw = nvpxcnw = {A s /E F )(T/E F ) K f- 2 (9) 

(dashed lines in Figs. 1131 and l*Hf|l does not exhibit such 
a maximum, but instead decreases monotonically with 
increasing doping. In Figs. Upland im we also plot A s /wo 
(thin solid lines), which shows that at low dopings, xss 
increases with increasing doping, despite the fact that 
the superconducting pairing strength A s decreases. 

We now consider the doping dependence of T c and the 
isotope effect on T c for a quasi- ID system that consists of 
an array of weakly coupled quasi- ID chains. We assume 
that the chains are spatially or dynamically fluctuating so 
that CDW order is dephased. The interchain Josephson 
coupling J is treated on a mean- field level^ so that T c 
is determined by the temperature at which 

2J XSS = 1 (at T — T c ), (10) 

where the numerical prefactor 2 is determined by the 
number of nearest-neighbor chains. Treating the inter- 
chain coupling with pcrturbative RG gives an equivalent 
result. Assuming J is doping independent, T c then ex- 
hibits a maximum at the same 5 where xss has a maxi- 
mum. 

The isotope effect exponent ax c is plotted versus <5 in 
Fig. El It is shown for various values of J = J/irv F , 
at fixed A and U . Unlike in BCS theory, ar c is not 
universal but depends on the interaction strengths and 
band-filling. However, qualitatively, it appears that the 
doping-dcpcndent behavior in which ar c is large near 
half-filling but decreases rapidly with increasing doping 
is generic (independent of interaction strengths and el- 
ph model). Note that if the parameters are tuned just 
right, ut c vanishes. A small ax c can even occur at the 
doping for which T c is maximum. Therefore, one should 
be careful not to assume that phonons are unimportant 
in unconventional superconductors for which ar c *C 1, 
such as in the cuprates at optimal doping. Figure El 
also shows oa,, which is weakly doping-dependent and 
negative. 

V. COMPARISON WITH OTHER WORK 

Below, we compare our results for the phase diagrams 
to some phase diagrams which have been previously com- 
puted. 

A. Half-filled extended Peierls-Hubbard model 

We first compare our results for the half-filled extended 
Pei-Hub model to recent Monte Carlo work by Sengupta, 
Sandvik, and Campbell^ on the same model. Following 
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FIG. 13: Doping dependence of the singlet superconducting 
susceptibility xss (thick solid lines) at T/loo = 0.1, CDW 
susceptibility xcdw (dashed lines) at T/uo = 0.1, and spin 
gap (thin solid lines), for the Pei-Hub model with U = 0.4, 
V = 0, Ef/loq = 5, and various values of Ap c i (labeled in 
plot). 
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FIG. 14: Doping dependence of xss (thick solid lines) at 
T /loo = 0.1, xcdw (dashed lines) at T /loo = 0.1, and A 3 /loo 
(thin solid lines), for the Hoi- Hub model [(a) and (b)], and the 
Pei-Hub model [(c) and (d)]. For all plots, U = 0.1, V = 0, 
and Ef/loo = 5. The values of Ah i and Ap e i are labeled above 
each plot. 



Ref. I22I we show a phase diagram in the 7-^0 A plane in 
Fig. 1161 This figure compares our result for the critical 
line for U/t = 2.5 and V/t = 0.625 (solid line) to the re- 
sult in Ref. [z^. (In order to plot our result in these units, 
we took the high-energy cutoff Ep in the RG theory to 
be t.) 

The quantitative disagreement near u>o/t ~ 1 can prob- 
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FIG. 15: Dependence of the isotope effect exponents qt c 
and oa, on the doping parameter S and interchain coupling 
strength J for the Hoi-Hub model with Ah i = 0.275 (a) and 
Pei-Hub model with A Pci = 0.25 (b). For both plots, U = 0.1, 
V = 0, and Ef/uq = 5. qa s is independent of J. 
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FIG. 16: Phase diagram of the half-filled Pei-Hub model, 
comparing of our result for the phase boundary (solid line) 
to the quantum Monte Carlo result of Ref. |22J (squares), for 
U/t — 2.5 and V/t = 0.625. The dashed line is our result for 
the same U/t but with V = 0. Here, 7 = -k (uo ~]i ) Ap e i /4 is 
the el-ph coupling constant in Eq. Q. 



ably be attributed to the fact that the assumption in the 
RG theory of a linear electronic dispersion becomes prob- 
lematic when ujq ~ Ep . At lower too/t, the agreement is 
excellent, especially considering the moderately strong 
value U = U/ (27rf) i=s 0.4. This gives one reason to be- 
lieve that the multistep RG method is, at the very least, 
qualitatively accurate for physically interesting values of 
[7-1. 

In Fig. we have also plotted our result for V — 
(dashed line) , which can be obtained from the simple an- 
alytic expression in Eq. I|59|l . The solid line in this figure 
is the only phase boundary in the present paper that re- 
quired numerical integration of the RG flow equations 
[Eqs. (|20|l and l|21|l ]. The phase boundaries in all other 
plots are given by analytic expressions derived in Section 

eh 




FIG. 17: Phase diagram of the spinless half-filled Holstein 
model from various authors, in the absence of el-el interac- 
tions. The solid line is an analytical result from Ref. l33t The 
line with circles and the line with squares denote the results 
of density matrix renormalization group (Ref. I20I) and two- 
step RG (Ref. list) , respectively. The line with triangles is 
the result from an exact-diagonalization method (Ref. l2ll) . 
(After Ref. H) 



B. Holstein model 

The phase diagram for perhaps the most interesting 
model studied in the present paper, the spinful incom- 
mensurate extended Holstein-Hubbard model, has not 
extensively studied in prior works. A much simpler re- 
lated model, which has been thoroughly explored, is the 
spinless half-filled Holstein model (without el-el interac- 
tions). We show the phase diagram for this model in 
Fig. El computed by various authors with a wide range 
of methods. Note that the result from the two-step RG 
technique (line with squares)^ is in good agreement with 
exact numerical methods. 

In order to see how Fig. 1171 changes when spin is in- 
cluded, we have computed a phase diagram in similar 
units for the spinful half-filled Holstein model in Fig. 1181 
including a Hubbard interaction. In the absence of el- 
el interactions, at weak el-ph coupling, the ground state 
of this model is always a CDW, for any finite u>q. This 
also holds in the strong coupling limit (AroI >• 1)£ In 
contrast, for the spinless half- filled Holstein model with- 
out el-el interactions, the transition to a CDW occurs 
at a nonzero value of <?, as shown in Fig. 1171 for weak 
coupling and proven in Ref. |9| for strong coupling. 

To study how Fig. ^| changes when the system is 
doped into the incommensurate limit, we have presented 
a diagram in the same units for the spinful incommen- 
surate extended Holstein-Hubbard model in Fig. [SJ In 
Figs. 151 and ITSl for technical reasons, we have defined 
g' = gy/Ep/vF, t' = t{E F /2t), and U' = U(E F /v F ). 

VI. METHODS AND DERIVATIONS 

Here, we provide more discussion of the technique and 
study the RG flows of the coupling parameters. We also 
derive the explicit expressions that were used to plot the 
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FIG. 18: Phase diagram in the g' /uJo-t' /loq plane of the spin- 
ful half-filled Holstein-Hubbard model for V — 0. The phase 
boundary is shown for various values of U'/uo, as labeled. 
The system is charge-gapped everywhere in the phase dia- 
gram, and spin-gapped in the CDW phase. For U' /too = 
(not shown), the ground state is a CDW for infinitesimal 
g'/wo. Here, g' » g, t' « t, and U' ~ U (see the text for 
their precise definitions). 



phase boundaries in the figures. 



A. Field theory for the 1DEG coupled to phonons 

To focus on the low-energy, long-wavelength physics, 
we work with continuum versions of Eqs. and 
The purely electronic part of our Hamiltonian is the stan- 
dard continuum model of the interacting 1DEG, in which 
the spectrum is linearized around the left and right Fermi 
points. The destruction field for fermions of spin a is 
written as a sum of slowly varying right and left moving 
fields: = e* kFX ip ha + e- ikFX ip- ha . The Hamiltonian 
density for the 1DEG, in the absence of el-ph interac- 
tions, is written as 7i = Ho + Wel-eb where the kinetic 
energy density is 



H = -iv F ^2 V^l,<ydxlp, h a, (11) 

r/,(7- ±1 

and the important short range el-el interaction terms are 
H e i_ e i = gi ^2 Ti^W-i^i^'^-i,^ 

+ 92 ^2 ^l,</^-l,<r'V ; -l,<T'^l,<r 

(T,<T'=±1 

+ 53 [e^-^* ^li^U-rV'i.-iV'i,! + H.c] 

+ 94 -o^-ffV^a- (12) 

We have assumed the system is spin-rotation invariant. 
The (72 and 54 terms describe forward scattering, the for- 
mer containing scattering on both left and right moving 
branches, and the latter containing scattering on only 
one branch. The g± term contains backscattering from 
one branch to the other. The 53 term contains Umklapp 



processes and is only important when 4fcp equals a recip- 
rocal lattice vector G; i.e., at half-filling (4kp — 2ir). For 
the extended Hubbard model in the continuum limit, the 
bare (unrenormalizcd) values of the g^s are given by 



5? 
9t 
9l 



- U-2V', 



5 3 

U + 2V', 
U/2 + 2V' 



(13) 
(14) 
(15) 



where the superscript indicates bare couplings and 
again V = — \^cos(2fcp), which equals V at half-filling. 

We incorporate el-ph interactions by defining retarded 
interactions 3i iPh , 52, P h, 53, P h, and g 4:Ph which play the 
same role as the g^s except that the energy transfer is 
restricted to be less than a cutoff uj C) which is approxi- 
mately the phonon frequency ojq (lu c will be defined more 
precisely below) . This corresponds to approximating the 
phonon propagator as a step function of frequency, which 
is a good approximation for the momentum-independent 
phonon dispersions we consider in this paper. In the Pei- 
Hub model, the bare el-ph couplings are given by 



53,pt 



51, P h - 

52, ph — 54, ph 



= -TTVpXpe 

= 0. 



For the Hoi-Hub model they arc 



5i, P h 



52, ph 



53,ph 



54, ph 



Hoi, 



(16) 
(17) 

(18) 



where Ap e i and AhoI are the positive, dimensionless cou- 
pling constants defined in Eq. (JSJ. Note that the cou- 
plings g? ph , g% ph , and g% ph are negative, indicating 
the attractive interaction induced by phonons. In the 
absence of el-el interactions, the sign of 3° ph is arbi- 
trary; likewise, the sign of g® is arbitrary in the absence 
of el-ph interactions. However, for the extended Hub- 
bard model coupled to phonons, once the sign convention 
03 = U — 2V' is chosen, it is required that g® ph > for 
the Pei-Hub model and 53 ph < for the Hoi-Hub model. 



B. Phase diagram of the 1DEG without phonons 

Before deriving the phase diagram including phonons, 
we briefly review the known quantum phase diagram of 
the 1DEG without el-ph coupling) 29 i 34 i 35 by identifying 
the conditions for various types of order to have divergent 
susceptibilities in the low-temperature limit. 

The sign of g± determines the existence or nonexistence 
of a spin gap: a 1DEG without el-ph coupling contains 
a gap to spin excitations for g^ < 0, and no such gap for 
5i > 0, regardless of band filling. A charge gap is only 
possible at commensurate fillings. 

An incommensurate 1DEG without a spin gap (Lut- 
tinger liquid), has a divergent 2k f SDW susceptibility 
when the Luttinger charge exponent 



K r = 



l27rvF + 2g2 + g 1 -2g° 2 



2nvi 



2g\ - g\ + 2g° 2 



(19) 
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TABLE I: Conditions for various ordering fluctuations to 
dominate for an incommensurate IDEG. See the text for the 
subdominant fluctuations. 



Type of order 



Dominates for" 



4/cf charge density wave (4k f) 
2/cf spin density wave 
2k f charge density wave 
Triplet superconductivity 
Singlet superconductivity 



A s = 0, K c < 1/3 
A s = 0, 1/3 < K c < 1 
A s > 0, K c < 1 
A s = 0, K c > 1 
A s > 0, K c > 1 



"For the Hubbard model, K c < 1 corresponds to repulsive inter- 
actions, K c > 1 to attractive interactions. The effect of a forward 
scattering el-ph interaction is to raise K c , while a backscattering 
el-ph interaction lowers K c and, if strong enough compared to the 
el-el repulsion, causes a gap in the spin sector (A s > 0). 



TABLE II: Conditions for SDW and CDW order in a half- 
filled, charge-gapped IDEG. 



Type of order 


Present for 


2k f spin density Wave 


A s =0 


2k f charge Density Wave 


A s > 



is less than 1, along with a logarithmically more weakly 
divergent 2k F CDW correlation. If K c < 1/2, the 4kp 
CDW correlation is also divergent; it is less divergent 
than SDW and 2k F CDW for 1/3 < K c < 1/2, but be- 
comes the dominant order for K c < 1/3. For K c > 1, 
a LL has a divergent triplet superconducting (TS) cor- 
relation and a logarithmically weaker divergent singlet 
superconductivity. 

An incommensurate spin-gapped IDEG (Luther- 
Emery liquid), is dominated by either SS or 2kp CDW 
correlations, depending on the value of K c . A LEL with 
K c < 1 is dominated by 2k F CDW. If K c < 1/2, 4k F 
CDW is subdominant, while if 1/2 < K c < 1, SS fluctu- 
ations replace 4k F CDW as the subdominant order. For 
1 < K c < 2, SS is the most divergent channel, and CDW 
is subdominant. For K c > 2, the only divergent corre- 
lations is SS. We summarize the phase diagram of the 
incommensurate IDEG in Table [I] 

At half-filling, a IDEG is charge-gapped for \g®\ > g® — 
2<?2, which for repulsive interactions is always satisfied. 
For g° < 0, the ground state of a charge-gapped system 
is an ordered, spin-gapped CDW (the Peierls instability) , 
otherwise, it is a SDW with no spin gap. 



C. The multistep RG technique 

To determine the phase diagrams, we apply the known 
one-loop RG flow equations^ 



dg 



2 dg c 

dl ~ 9l ' dl 



2 d 93 

~ 53 ' ~df 



dg±, P h 
dl 

rfff2,ph 

dl 



''.I ,ph \^9l ± '/:'. 

rfff4, P h _ dm _ „ 

dl dl ~ ' 



1 



doj c . , dfi 

— =u c ( %VF+g+ , ph ), _ 



j.9c + g±, P hJ ,(21) 
(22) 

ttvp/j, (23) 



where we defined 



9c = 9i - 2g 2 , 

.9±,ph = 91, ph ± 33,ph, 

I = (■kv f )~ 1 \ti(E f /uj), 



(24) 
(25) 
(26) 



and u) is the running cutoff. The above expressions apply 
for E F lo ^> u>f), /i. If, instead, ^ 3> u> lu q , the same 
equations apply, but with 53 = <73 iP h = 0. ^From Eq. 
(|20|l . we see that a repulsive gi renormalizes in the same 
way as the Coulomb pseudopotential in a Fermi liquid-it 
is scaled to smaller values as one integrates out high- 
energy degrees of freedom. However, in ID the backscat- 
tering and Umklapp el-ph interactions <?i, p h and fl^ph 
are strongly renormalizcd. Furthermore, there are cross 
terms <7±. p h<7i in Eq. I|21|) . which means that the RG 
flows of <7i, p h and g^.ph are strongly influenced by direct 
el-el interactions. 

The two-step RG procedure is as follows: Assuming 
A s < ojq and the system is at half-filling (p = 0), we first 
integrate out fermionic degrees of freedom between the 
high-energy scale Ep and the phonon energy uiq using 
Eqs. (|2*0"|l and (f2~T|) . Once ujq is reached, total effective 
interactions g* ot (wo) at this energy scale are determined 
by adding the effective el-el coupling to the effective el-ph 
coupling 



9i Ot (u ) = 3i(w ) + &, P h (<*>(>)• 



(27) 



~939c 



(20) 



Below this energy scale, there is no difference between re- 
tarded and instantaneous interactions, so g\ ot renormal- 
izes as a non-retarded interaction using Eq. 1)20(1 with 
gl ot (too) as the initial value. If the system is far enough 
away from half- filling such that /1 ~ Ep, the method 
is identical except that one sets g§ = g® ph = at the 
start. For the above cases, it is clear that the renormal- 
ized couplings <7' ot (a;o) and the renormalized cutoff loq 
play the same role in the ID electron-phonon system as 
<7° and Ep, respectively, do in the pure IDEG. Therefore, 
to determine the phase diagram of the IDEG coupled to 
phonons, we can use the known phase diagram of the pure 
IDEG, and simply replace the <?°'s there with <7* ot (wo)'s. 
For more general fillings (0 < fi < Ep), a three-step RG 
method is necessary, which we will elaborate at the end 
of this section. 

Note that luq is the physical phonon frequency, which 
is related to oj c by the expression ujq = lu c (luq), where 
lo c (loq) is the renormalized value of u) c at the energy scale 
ujq, determined by Eq. (|23|l . Likewise, we define fi as 
the physical value of the chemical potential relative to 
its value at half-filling; the bare chemical potential fi° is 
chosen such that it flows to the value /j, after integrating 
out degrees of freedom between Ep and fi. 
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D. Incommensurate filling 

We first derive the scaling of the coupling constants 
and compute the phase diagrams for the case when the 
system is doped far into the incommensurate limit (/i ~ 
E F ). 



1. Scaling of the coupling constants 

Using Eqs. and J2IJ with g 3 = g 3iP h = 0, we 

integrate out degrees of freedom between Ep and u> with 
Ep > oj > u>o, to obtain the follow effective couplings for 
incommensurate systems: 



5eM 
91 (w) 

5i, P h(w) 
where 



.9c 



9°i 



1 + tftV 




gi(u) \ 3/2 ( E F 



1 exp(-gg a ;/2) 

(1 + g \ x )m 



(28) 
(29) 

-g°/2nvF 

J 

(30) 
(31) 



and again g° = gi(E F ). 

In the absence of el-el interactions, the el-ph backscat- 
tering coupling gi p h flows to stronger values according 

to ffi,phH/5i, P h = (1 - l^phlO" 1 - In Fig. d the im- 
portant influence of extended Hubbard interactions on 
5i,ph(w) in studied. In Figs. I19f a) and I19f c1. we show 
the scaling for several fixed finite values of U/\g® ph |. Fig. 
I19f e) shows the scaling in the limit U 3> \g® ph \, which is 
given by 

Si,phM/3?, ph w HUl), (U » \gl ph |), (32) 
where for future convenience we define 

exp(G c a;/2) 



h(x) 

d 
G, 



(1 + Gix) 3 / 2 ' 
l-2(V/U), 
1 + 6(V/U). 



(33) 

(34) 
(35) 



^From Fig. ^3 we see that the flow of <7i, pn is very sen- 
sitive to the parameter V/U. If V/U < 1/6 and U/\g° lph \ 
is large enough, then the g±, P hffi term in Eq. (|2"T|) causes 
ffi,ph(w) to initially flow to weaker values. However, if 
V/U > 1/6, then gu P h(uj) / g^ ph > 1 for all uj < E F , re- 
gardless of U/\g° h |. The driving force for this increase 
in |ffi, p h(w)| is the 5±, p hffc term in Eq. (gU- 

To better understand why increasing U can sometimes 
enhance superconductivity at small U (see Fig. |SJ), we 
study the dependence of <7i jP h(wo) on U at fixed Ef/loq 
in Fig. 1201 In this plot, for small U, the effective t/i, p h is 



Incommensurate filling 

(Any model): 



u/ = - 5 



Half-filling 

(Model as shown): 

u/|g °| = o.s 



2 - 







' (a) 




' 7 / 






■ / s 
























V/U 


= 




V/U 


= 1/10 




V/U 

1 


= 1/6 



2 - 



s 1 



(b) 



i / 
i I 



Half-filled Pei-Hub 

(V = 0) 
Half-filled Hoi-Hub 

(V = 0) 

i I i 



0.25 
Ul 

U/ lO = 2 



0.5 



0.25 
Ul 



0.5 



2 - 




y/u = o 

— y/G = 1/10 

V/U = 1/6 



"T- 



u/|g,, PB °[ = 2 



(d) 



/ 
/ 

/ 

/ Half-filled Pei-Hub 

(V = D) 

Half-filled Hoi-Hub 

(V = 0) 



1 

Ul 



u>> iO 



1 

Ul 



2 - 



-V/U = 

wu = 1/10 

-V/U = 1/6 



2 - 



(f) 



Half-filled Pei-Hub 

(V = 0) 
Half-filled Hoi-Hub 

(V = 0) 



1 

Ul 



1 

Ul 



FIG. 19: Dependence of the effective backscattering el-ph 
coupling (?i !P h(uj) on Ul = (U/ttvf) ln(-Ep/a/) for u > wo- 
Plots (a), (c), and (e) are for systems doped into the incom- 
mensurate limit, while (b), (d) and (f) are for the half- filled 
Pei-Hub and Hoi-Hub models with V — 0. For each plot, the 
ratio of U to the bare el-ph coupling g° ph is held fixed at the 
value indicated above the plot. For each curve in (a), (c), and 
(e), the ratio V/U is held fixed at the value indicated. 



0.8 
>"-0.6 
£0.4 
? 0.2 

0.0 



Incommensurate filling 
X = 0.2, E F / M[1 = 5 



-V/U = 

y/u = 1/10 
-y/y = 1/6 

-V/U = 1/4 



(a) 




0.0 0.1 0.2_0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 
U U 

FIG. 20: Dependence of gi, p h(wo) on U for incommensurate 
systems with (a) Ef/ujo = 5 and (b) Ef/uJo ~ 10 (b). The 
horizontal dotted line indicates the value of the bare el-ph 
coupling A = — glph/nVF = 0.2, where A stands for either 
AhoI or Ap c i, depending on the model. For each curve, the 
ratio V/U is held fixed at the value indicated in (a). 
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stronger that its bare value (<7i,ph( w o)/<?i ph > 1)- How- 
ever, for V/U < 1/6, increasing U causes |pi, p h(wo)| to 
decrease. This can cause Kf* and xss to increase, as 
in Fig. [S] However, if we set V/U > 1/6 and hold this 
ratio fixed while increasing U, then |5i, p h(^o) I increases, 
causing K f 1 and xss to be suppressed. 

Note that in real materials, one typically expects 
V/U > 1/6, in which case | <7i, P h(<-^o) | > |ffi p hl or even 
lffi,ph(w )| > |ff° ph l- Tne requirement for a spin gap 
is |fl , i, P h(wo)| > 5i(wo)j which can be achieved, in many 
cases, with even a small amount of retardation. (Note 
that for repulsive el-el interactions, < .gi(^o) < <7i-) 
Therefore, it is possible for a slightly retarded el-ph in- 
teraction to create a divergent superconducting suscepti- 
bility even when the bare interactions are predominantly 
repulsive (|ff? jph | < #?). 

2. Luttinger liquid to Luther-Emery liquid transition 

The phase boundary between the LL and LEL phases 
is given by the condition g\ ot (uJo) — 0. This condition 
determines the following critical value of the bare el-ph 
coupling 



\ Gap 
A Pei 



A 



Gap 
Hoi 



u 



exp{G c Ul /2) 
Giy/l + GiUl 



f(Ulo) 



(36) 



(LL— LEL transition), 



where we define 



For Apoi > A p ^ p or Ah i > ^Hof> tne ground state of 
the incommensurate 1DEG is a spin-gapped LEL; other- 
wise, it is a gapless LL. The condition A = A p ^j p = A^ p 
therefore determines the transition lines in Figs. H an d 
121 as well as the thick solid lines in Figs. |3 01 and 
Note that in Fig. |3 we used that the "scaled" critical 
couplings A^j P Zo and X p ^Iq depend only on the two pa- 
rameters UIq and V/U. 

3. Susceptibilities and spin gap in the LEL phase 

In the LEL phase, the potentially strongly divergent 
part of the low-temperature susceptibilities for SS and 
2kp CDW are given by Eqs. I© and ©, respectively, 
where A° ff is the effective Luttinger charge exponent af- 
ter integrating out states between Ep and uiq, given by 



f(y) 



dx h(x) 



(37) 



A" 



off 



l2TTV F + 2gl ot + g t ° t {u ) 



2itvf 



i\ ot 

2 5 | ot 



where 



ff * ot M = g c + .gi,p h M-2 32 ° 

tot _ , 

9\ — 9i t .94 iP h- 



. ph > 



(38) 

(39) 
(40) 



Integrating out states below the energy loq does not fur- 
ther renormalize g* ot (and therefore A c ). Note that the 
effective charge and spin velocities are also renormalized 
due to phonons, the former of which is 



(27T)" 



[2ttv f + 2g\ ot ] 2 



(41) 



For —1 < g\ ot (ujo) < 0, A s is given approximately 
by the energy scale below ujq at which the RG analysis 
breaks down because the effective g\ ot /TTVF has grown to 
— 1. Since 



<?i ot M = 

this gives 



ffi ot M 



l + [g\ ot {uj )/T:v F ) In(w /w) 



for lo < loq , 

(42) 



a; eexp[-7r^/|^ ot (w )|]. (43) 
as u> — > 0; therefore, A s = 



For g\ ot (u) ) > 0, g[ ot (uj) 
0. 



4- Competition between SS and CDW in the Hoi-Hub model 

The thin solid line in the incommensurate extended 
Hoi-Hub model phase diagrams of Figs. |3 0J and 03 
which we call the "superconducting transition," is de- 
termined by A° ff = 1. This condition is satisfied for 



AhoI = Atr 3 .'," 1 " and AhoI — A 



ss. 



where 



,SS,± _ j-j 
A Hol - U 



H±JH 2 



G c 



2f(Ul ) 



(44) 



(superconducting transition) 



and 



H 



h(Ul ) , G c 



4/([// ) 



(45) 



Therefore, as long as the square root is not imaginary, 
for fixed Ulo > 0, there are two critical values of the 
bare el-ph coupling determining the boundary between 
the phase with dominant CDW and the phase with dom- 
inant SS. The most divergent correlation is SS provided 

Note that 



G c < 2f{Ul )H 2 and A| s ol ~ < A Ho i < A 1 
divergent a TS susceptibility is not present in any phase 
diagrams, since for repulsive el-el interactions, A° ff < 1 
whenever g\ ot (u!o) > 0. 

The phase boundary shown as a dashed line in Figs. 
01 EI and is defined by Kf = 1/2. For the extended 
Hoi-Hub model, this condition occurs at AhoI = A Hoi W 
with 



v ss,+ 

Hol 



iCDW 
^Hol 



= U 




(46) 



(Kf = 1/2 transition), 
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where 



A-5h(Ul Q ) Q 
' = 8/(m ) 4' 

= 3/u-9(y/u) - 1. 



(47) 
(48) 



For AhoI > -^Hoi W ' a divergent SS susceptibility is not 
possible. 



5. Competition between SS and CDW in the Pei-Hub model 

For the extended Pei-Hub model with repulsive el-el 
couplings, K° s < 1 always, which means a phase with 
dominant SS is impossible. For this model the condition 
Kf = 1/2 occurs at the critical el-ph coupling value 



iCDW _ fr 
A Pei — u 



5h(Ul ) 
2Q 



+ f(Ul ) 



(49) 



(Kf = 1/2 transition). 



For A Pci > Ap^ w , the SS susceptibility is never divergent, 
on 

in Figs. EUc) and Eld). 



The condition Ap i = Ap^ w determines the dashed line 



E. Half-filling 

At half-filling, Eqs. l(2T?|) and (J2TJ can be integrated 
analytically if one takes gj? = —g®. For that case, 



3+,phM 

where 
X -- 



-93 (w) 



9" 



l + g°J 

g+, P h \ ( giH 



5 +,ph 



3-,ph 




3? 



5cM 



(50) 



-1/2 



, 3? 



3/2 



(51) 

/ \\ 3/2 



dx{l+glx)-^(l + glxyi\ 



Y= / ^(l+ ff ? a; )- 3 / 2 (l+ 5 °a ; )- 3 / 2 
Jo 

2 [ o , o 3? +9° C +H9° c l 
3i t 3 C 



(3? - 3c°) 2 



V(l+3i0(l + 3 c °0 



(52) 



(53) 



(54) 



and gi(u>) is given again by Eq. (|29(l . Since our assump- 
tion of g° = — <?3 is satisfied for the Hubbard model (with 
V = 0), the above results determine the scalings of <7i jP h 



and 33, p h for the Hoi-Hub and Pei-Hub models, which 
are 



giM = 53 < phM = (i+^j V (TTW 

(half— filled Holstcin— Hubbard model), (55) 



ffx lP h(w) = ~33,ph(w) 



3l,ph 



(half— filled Peierls— Hubbard model), (56) 

where X' and Y' are the values of 2X and 2Y, respec- 
tively, for the case = — 3c = ^ : 



A " = Z7 



2-2 



2/ 



1 - t/7 
l + [/Z 



arcsin(J77) 



(57) 
(58) 



In the absence of el-el interactions, for either model, 
3i,ph(w) increases in strength as I is increased according 
to 3i,ph(w)/3i lPh = (! - %i )P hlO -1 - As shown in Fig. 
IToT bl. (d), and (f), turning on a repulsive U has the op- 
posite effect for the half-filled Pei-Hub model compared 
to the half-filled Hoi-Hub model: for the former gi jP h(u;) 
increases even more rapidly with increasing I than be- 
fore, while for the later (7i. p h(w) increases less rapidly 
than before. This is due to the bond (site) centered na- 
ture of the Peierls (Holstein) el-ph interaction, and shows 
up as a sign difference in g® ph for the two models. 

For U > |#° ph |, the behavior of g hph (u)/g^ ph is given 
by the square roots in Eqs. (|5*5|l and (55} , an( i is shown in 
Fig. I19f f). In this limit, for the Hoi-Hub model, gi iP h(w) 
flows to weaker values. Since the spin gap is enhanced 
by a strong <?i. p h, and the charge gap is enhanced by a 
strong 53 iP h, we see that the off-diagonal phonon mecha- 
nism in the Pei-Hub model is more effective in enhancing 
both the charge and spin gap compared to the diagonal 
mechanism in the Hoi-Hub model. 

At half-filling, the transition line between CDW and 
SDW phases, called the Mott-Peierls transition, is deter- 
mined by g\ ot (uJo) — 0. The critical el-ph couplings that 
define this phase boundary are then 



\ GAP 
A Pci 



Uy/1 - {Ulp) 2 

(1 - Ul )- 1 + 2Ul 



(59) 



(Peierls— Hubbard model: Mott— Peierls transition), 



xGAP _ fj 
A Hol — u 



4~3W ! ^° - 2 arcsin(W ) 
V i + u <o 



(60) 



(Holstein— Hubbard model: Mott— Peierls transition). 
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For Ap c i > Ap^ p in the half-filled Pei-Hub model, or 
AhoI > ^h^ P m the half-filled Hoi-Hub model, the ground 
state is a spin-gapped, ordered CDW. We plot the tran- 
sition line given by A Pci = A^ p in Figs. H2» and EH 
and the line determined by AroI = ^h^ P hi Figs. IT2T a') 
anddl 



F. Near half-filling 

Using the two-step RG technique, we have derived 
phase boundaries for the strongly incommensurate case 
/i ~ Ep, as well as the half-filled case \i = 0. For the 
more general case < fi < Ep, in other words at fill- 
ing near but not equal to half-filling, a three-step RG 
method is necessary. The three distinct crossover scales 
are the high-energy Ep, low-energy scale ujq, and chem- 
ical potential /i. As before, retarded interactions only 
renormalize when integrating out states between Ep and 
luq. However, now, 173 and g3 iP h only play a role when 
integrating out states at higher energies than \x (and if 
fi < luq, for states between /j, and ujq, only plays a 
role). 



for either model, where 
r sui 

Z = 



(67) 



Since g\ renormalizes in the same way for the half-filled 
and incommensurate cases, we can just integrate from 
Ep to u in one step using Eq. (J23). 

Again, the condition g\ ot (uJo) = 9x(^o) + 5i, P h(w ) = 
determines the transition to a spin gap, which leads to 
the critical values 



Pci 



gyg - (cWo) 2 ] 3 
4 + £ + 2cUl [3- (cUl ) 2 } 

(Peierls— Hubbard model), 



(68) 



\gap 
A Hol 



u 




1 - cUl 
{1 + cUloY 



2 arcsin(cC// ) 



(69) 



(Holstein— Hubbard model), 
where we defined 



1. Doping dependence of the phase boundaries 

We now employ the three-step RG technique to derive 
the doping dependence of the phase boundaries for V — 
0. First consider the case ujo < H < Ep. We begin 
by integrating out degrees of freedom between Ep and 
/i = uJo(Ep/ujo) s , resulting in an effective 31, p h of 



ttvfX 



Hoi 



1 - cUl 



l-XaaX/Uj V {1 + cUlof 



TTOpApei \ 
1 - 



(61) 



(62) 



for the Hoi-Hub and Pei-Hub models, respectively, with 
c = 1-5, (63) 



X 



Y = 



'2 2\l - — — arcsin(c[/^ ) 



1 + cUlo 



2cUh 



y/1 - (CUlo) 2 



(64) 
(65) 



Next, <7i lP h(A < ) is used as the initial value to integrate 
from fi to ujq, employing the RG flow equations without 
g 3 and fl3, P h, resulting in 



ffl, P h(wo) 



gi.phQ") ) / cxp(8Ul ) 



l + gi . ph (fi)Z/Uj V (l + <^o) 3 



(66) 



£ = (l + Ul )F-4{l + cUl )+Z, (70) 
F = e sm ° /2 (l + SUl )- 3/2 . (71) 

The system is spin-gapped for Ap c i > Af,^? or Ah i > A§^ 
for the Pei-Hub and Hoi-Hub models, respectively. The 
phase boundaries determined by Ap C i = Ap^ 5 and Ah i — 
A™ are shown in Figs. |Hl-^2 as thick solid lines. For 
8 — 1, App? = A^, and we recover the incommensurate 
LL-LEL transition [Eq. wit h V = 0. 

We obtain analytic expressions for the remaining phase 
boundaries by requiring that K° s equals 1 or 1/2 (de- 
pending on the phase boundary), using 



Vo) = 9c(p) + £fi, P h(wo) - 2gl 



P h 



(72) 



with (? c (/i) = — £7/(1 — cUIq). The results for the critical 
couplings, for ujq < fi < Ep, are then 



A Hol - U 



\cdw j r 

A Hol — U 



B ± y/B 2 - SAC/2 



\ cdw 
A Pci 



u 



D + ^D 2 + 5SAE/4 

Y F/E + Z 

- (Mo) 2 ] 3 



-1 -1 



(73) 
(74) 

(75) 



with the definitions 



A = 


(sx + zy\ 




(76) 


B = 


[(2S-F)A + C]/4, 




(77) 


C = 


{l-cUh)-\ 




(78) 


D = 


A[S{A~hEX) - 5(EZ 




(79) 


E = 


(6/[7 + 3)/5-C, 




(80) 


S = 


(l + cUl ) 3/2 (l-cUl ) 




(81) 
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For 5 = 1, Eqs. (J73J), O, and J7SJ) reduce to Eqs. (gU, 
(|4fc>p . and l|49l) . respectively, with V = 0. The conditions 

AhoI = ^Haf ano - AhoI = ^HoT determine the thin solid 
line in Figs. l^land irTT aL In Fig. Illf bl. U is large enough 
such that SAC > 2B 2 everywhere in the plot, therefore 
the phase with dominant SS is not present. The condition 
AhoI = Ag^Y determines the dashed line in Figs. ISI HTT aL 
andEtb). The dashed line in Figs. UHI Etc), andEtd) 
are determined by Ap i = Apgf. 

For the case /i < ujq < Ep, since g\ scales in the same 
way at all energies and does not depend on g^, the spin- 
gap phase boundary is independent of /i, and is given by 
Eqs. (HE) and with S = 0, or by Eqs. ^ and ffi . 



2. Doping dependence of susceptibilities and isotope effects 

The doping dependence of the susceptibilities in the 
LEL phase (Figs. H3land lT^jl is computed with the three- 
step RG method using Eqs. iJBJ and 10, combined with 
Eqs. (S§, g3J, and J72J|. 

The isotope effect on T c (Fig. I15|) is computed via 



OtT a 



1 AT C /T C 
"2 AC/C 



(82) 



where C = E F /uj = e l ° and A( < (• Here, AT C = T c ' - 
T c , where T' c is the transition temperature determined 
from Eq. JT0|) after changing l — > l' a = \n(( + A£), 
c —> c' = cIq/Iq, and 6 — > 1 — c'. The changes in c 
and 5 are required so that, when changing £, only the 
energy scale u>q is changed, and the energy scales Ep and 
remain fixed. The isotope effect on A s is determined 
in a similar fashion. 



VII. CONCLUSIONS 

We have explored the influence of the el-ph interaction 
on the quantum phase diagram of the most theoretically 
well understood non-Fermi liquid, the interacting 1DEG. 
The backward and Umklapp scattering portions of the 
el-ph interaction are strongly renormalized, often toward 
stronger couplings. Even in the presence of strong el-el 
repulsion, a weak, retarded el-ph interaction is capable of 
creating a spin gap and causing divergent superconduct- 
ing and/or CDW susceptibilities (true long range order 
is formed when weak coupling between ID chains is in- 
cluded). The ground state is strongly dependent on the 
band-filling, and, especially at or near half-filling, depen- 
dent on the microscopic model of the el-ph interaction. 
Compared to higher dimensions, the zero-temperature 
phase diagram is far more complex, and, away from com- 
mensurate filling, contains a subtle competition between 
SDW, CDW, and superconductivity. The fact that direct 
el-el interactions strongly influence the renormalizations 
of the el-ph interactions adds to the richness of the phase 
diagram. In ID, intuitive concepts that apply to higher 



dimensional Fermi liquids, such as the suppression of su- 
perconductivity by repulsive interactions at weak cou- 
pling, must sometimes be abandoned. 

When the bare el-el repulsion is much stronger than 
the bare el-ph induced attraction, in ID, unlike in higher 
dimensions, it is not a requirement that Ep 3> u>q for 
the superconducting susceptibility to diverge. (In fact, 
in ID, very large values of Ep/luq are harmful to su- 
perconductivity.) Note that in the high-temperature su- 
perconductors, where Ep /coq ~ 5 and the el-el repulsion 
is strong, it has been correctly argued that the small 
value of Ep/uo rules out conventional phonon-mediated 
superconductivity^ It is interesting to point out that the 
arguments there apply only to a Fermi liquid and not the 
quasi-lDEG. 

We now qualitatively summarize the phase diagrams, 
for the case of repulsive el-el interactions, beginning with 
a system that is far from half-filling. In this case, the 
charge sector is gapless. For either the Hoi-Hub or Pei- 
Hub model, the spin-gapped LEL phase is favored by 
small U, large V, large A, and large Ep/iv . The LL 
phase is favored by large U, small V, small A, and small 
Ep/ujQ. For the Hoi-Hub model, a dominant supercon- 
ducting fluctuation is favored by small U, small V, mod- 
erate AhoI, and small Ep/luq. For either model, the phase 
with dominant 2kp CDW and subdominant supercon- 
ductivity is favored by moderate Ep/uo and moderate 
A (the dependence on U and V is more subtle-see the 
phase diagrams). For either model, the phase with dom- 
inant 2kp CDW and subdominant Akp CDW is favored 
by large V, large A, and large Ep/uJo (see the diagrams 
for the subtle dependence on U). 

Moving the incommensurate system toward half-filling 
increases the stability of the LEL phase relative to the LL 
phase in the Pei-Hub model, but decreases the stability 
of the LEL phase relative to the LL phase in the Hoi- Hub 
model. In the Hoi-Hub model, moving toward half-filling 
suppresses the phase with dominant superconductivity. 
For both models, moving toward half-filling decreases the 
stability of the phase with subdominant SS and increases 
the stability of the phase with subdominant Akp CDW, 
but this effect is more pronounced in the Pei-Hub model. 

In the Hoi-Hub model at half-filling, a spin-gapped 
CDW phase is favored by small U and large AhoI- The 
SDW phase with no spin gap is favored by large U and 
small AhoI- The half- filled Hoi- Hub phase diagram is 
weakly dependent on Ep/luq compared to the Pei-Hub 
phase diagram. In the half-filled Pei-Hub model, the 
spin-gapped CDW phase is favored by large Ap e j, large 
Ep/u)Q, large V, and any U other than moderate val- 
ues. The SDW phase is favored by small Ap e i, small 
Ep/uio, small V, and moderate U. Both models are 
charge-gapped at half-filling. 

We have studied the strong doping dependencies of the 
phonon-induced spin gap and various susceptibilities in 
the Luther-Emery liquid phase. The spin gap and charge 
density wave susceptibilities decrease monotonically as 
the system is doped away from half-filling. However, 
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the superconducting susceptibility, and therefore T c in 
a quasi-lD system with fluctuating chains, can vary non- 
monotonically with doping and exhibit a maximum at 
some "optimal" doping. 

Partially motivated by the unconventional doping- 
dependent isotope effects observed in the cuprate high- 
temperature superconductors, we have computed isotope 
effects in the quasi-lDEG coupled to phonons, since it is 
perhaps the most easily studied unconventional phonon- 
mediated superconductor. The calculated isotope effects 
bear a qualitative resemblance to those observed in the 



cuprates, as summarized in Section [I] 
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